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ABSTRACT 

We examine the rest-frame far-infrared emission from powerful radio sources with 
1.4GHz luminosity densities of 25^1og(£i.4/WHz~ 1 )^26.5 in the extragalactic Spitzer 
First Look Survey field. We combine Herschel/ 'SPIRE flux densities with Spitzer /IRAC 
and MIPS infrared data to obtain total (8 — 1000 /im) infrared luminosities for these 
radio sources. We separate our sources into a moderate, 0.4<z<0.9, and a high, 
1.2<z<3.0, redshift sub-sample and we use Spitzer observations of a z<0.1 3CRR 
sample as a local comparison. By comparison to numbers from the SKA Simulated 
Skies we find that our moderate redshift sample is complete and our high redshift 
sample is 14 per cent complete. We constrain the ranges of mean star formation rates 
(SFRs) to be 3.4-4.2, 18-41 and 80-581 M Q yr- 1 for the local, moderate and high 
redshift samples respectively. Hence, we observe an increase in the mean SFR with 
increasing redshift which we can parameterise as ~ (1 + z)® , where Q = 4.2 ± 0.8. 
However we observe no trends of mean SFR with radio luminosity within the mod- 
erate or high redshift bins. We estimate that radio-loud AGN in the high redshift 
sample contribute 0.1— 0.5 per cent to the total SFR density at that epoch. Hence, if 
all luminous starbursts host radio-loud AGN we infer a radio-loud phase duty cycle 
of 0.001-0.005. 

Key words: galaxies: active, star forming, infrared: galaxies, radio continuum: galax- 
ies. 



1 INTRODUCTION 

There is now strong evidence that powerful active galactic 
nuclei (AGN) played a key role in the evolution of galax- 
ies. T he correlation of centr al black hole and stellar bulge 
mass (Ma gorrian et al. 1998), and the increased prevalence 
of sta r formation ( Hopkins fe B cacom 2006; Giavalisco et al.l 

I l2010h 



2004) and AGN activity (|Wall et al.ll2005l ; lAird et al 



at earlier epochs suggest that the growth of the black hole 
is somehow related to the growth of the host galaxy. In 
the local Universe we see little evidence of high star forma- 
tion rates (S FRs) in galaxies with powerful radio-loud AGN 
activity (e.g ICondon et all 1 19981 : iMauch fc Sadler! l2007h . 
In the distant Unive rse, z > 1, luminous ra dio galaxies 
l|Sevmour~ et al. 200?]) and powerful starbursts (|Borvs et al.l 
l2005l : ICasev et alj2009h are both hosted by massive galaxies, 
suggesting a common parent population. The idea that these 
processes are likely to be connected at the epoch when black 
holes and galaxies went through their most rapid phases 
of growth has been invoked within various semi-analytical 
mod els in order to reconc ile those models with observations 
(e.g. ISpringel et ai1l2005h . 

This connection between central black hole growth and 
star formation rate is often considered in the context of 'feed- 
back' process(es), as the former is postulated to regulate the 
latter. In particular there is observational evidence, as well 
as theoretical models, in which the jet from an AGN can 
produce either positive or negative feedback, where the jet, 
traced by its radio emission, stimulates or quenches star for- 
mation respectively. There is some observational evidence 
of positive feedback, whereby star formation is triggered by 
an A GN jet, e.g. in Minkows k i's Object by a je t from NGC 
541 |van Breugel et all 1 19851 : ICroft et allhood ). as well as 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA, 
f E-mail: npsOOmssl .ucl . ac .uk 



theoretical models which suggest that the shocks generated 
by jet propagation can trig ger collapse of over-dense clouds 
and l ead to star formation (| Fragile et a l. 2004; S axton et al.l 
2005). Negative feedback by AGN jets would likely require 
the removal of fuel for star formation, evidence for which 
are the powerful AGN-induced outfl ows which have been 
seen in high redshift radio galaxies l|Nesvadba et al.l 120061 . 
2008). Such a scenario has also been proposed to regulate 
the growth of mas s ive galaxies in se mi-empirical models 
ICroton et al.ll2006l : iBower et af]|2006l ). but this process is 
only important globally at late times, z < 1. At earlier times 
it would be most important in halting the growth of the most 
massive galaxies. 

Star formation in powerful AGN has been difficult to 
trace so far. This difficulty is due to heavy contamination 
in traditional diagnostics by emission from the AGN (e.g. 
UV luminosity or optical emission line strengths) as well 
as obscuration by gas and dust. However the far-infrared 
(far-IR) presents a window in the electromagnetic spectrum 
where AGN emission is weak and star formation, if present, 
can dominate. AGN dust emission tends to peak in the 
near/mid-IR so far-IR emission should be a cleaner measure 
of SFR than other traditional methods. It is also possible 
to use the near/mid-IR to model and subtract any potential 
AGN contribution to the far-IR (e.g. lHatziminaoglou et al.l 
l2008h . 

There is evidence for extreme SFRs in many power- 
ful high redshift radio galaxies (z > 2, 1.4 GHz luminos - 
ity densities, L 1A > 10 27 WHz~\ lMilev fc De Breuck|[2008 



from their strong sub-mm emissi on ([Archibald et alfeoOll ; 
iReuland et al.ll2004 : [Greve et al.l 12006^ 7 their mid-IR spec- 
tra ( Seymour et al. 120081 . J. Rawlings. 2011, in prep.) and 
the spectacular (>100kpc) Lyq haloes someti mes observed 
j|Reuland et al.ll2003l : Fvillar-Martfn et al.ll2003t ) showing the 
extended gas that can provide the fuel for star formation. 
To compliment future targetted Herschel studies of the rare, 
very powerful radio- loud AGN, we examine in this work less 
luminous radio-loud AGN, 26.5 > log^Li^/WHz -1 ) > 25, 
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which can be found in reasonable abundance over areas of 
a few square degrees. We use this definition of 'radio-loud' 
AGN , based on radio luminosity density {e.g., iMiller et al.l 
1990), in order to avoid making any distinction between type 
1 and type 2 AGN, i.e. AGN classification based upon op- 
tical spectroscopy, where different amounts of AGN obscu- 
ration may affect the relative amount of optical emission. 
Although, as we shall show, most of these sources are also 
' radio -loud' when using the definition of iKellerman et al.l 
(1989, 5 GHz over B— band luminosity > 10). Star formation 
in these less luminous radio-loud AGN remains poorly stud- 
ied, as there has been no systematic follow-up of such sources 
above z > 0.1. Recently, the importance of rad io-loud AGN 
in thi s luminosity range was demonstrated by ISaiina et al.l 
(2007) who found that 40 per cent of z ~ 2 ULIRGs with 
deep silicate absorption features were radio-loud and those 
authors postulated that such sources are transition 'feed- 
back' objects after the radio jet has turned on, but before 
feedback has halted black hole accretion and star formation. 



The SPIRE instrument jGriffm et alj|20ich on board 
the Herschel Space Observatory ( Pilbratt et al.l l2010Tl gives 
us a clear view of the far-IR/sub-millimeter Universe at 
wavelengths where many galaxies emit most of their lumi- 
nosity. The Herschel Multi-tiered Extragalactic Survey (Her- 
MESJ, Oliver et al. 2011, in prep) provides deep infrared 
SPIRE data over many of the best studied extra-galactic 
survey fields. Recent results from Herschel show that SPIRE 
detected AGN in deep HerMES fields have far-IR colours 
similar to the bulk of the SPIRE popula tion which are be- 
lieved to be star formation dominated l|Elbaz et al.l |2010| ; 
lHatziminaoglou et al.l I2010T ) and modeling of their spec- 
tral energy distributions (SEDs) suggests the SPIRE emis- 
sion in AGN is dominated by a star forming component 
|Hatziminaoglou et alj|201oh . 

The work presented here uses Herschel /SPIRE observa- 
tions of the Spitzer Extragalactic First Look Survey (FLS) 
field taken as part of the Herschel Science Demonstration 
Phase (SDP) in October to November 2009. Of the fields 
observed in SDP this field had the best combination of wide 
area, uniform radio coverage and good multi-wavelength 
follow-up. We present our sample of moderate and high 
redshift radio-loud AGN in §2 and, we derive IR luminosi- 
ties and star formation rates in §3. We present our re- 
sults in §4 and discuss them in §5. We conclude this pa- 
per in §6. Throughout we use a 'concordance' cosmology of 
f2 M = 1 - f2 A = 0.3, Q = 1, and H = 70 kms" 1 Mpc" 1 . 



2 SAMPLE 

2.1 Radio Sample and Cross Identification 

Our radio d ata come from the 1 .4 GHz Very Large Array 
catalogue of lCondon et~ai1 |2003) which is complete down 
to 0.115 mjy (5c). We restrict our analysis to a region of 
the FLS with complete optical and near/mid-IR coverage, 
defined by 257.8° < RA < 261° and 58.6° < dec. < 60.4°. 
These optical to mid-IR ancillary data were taken from 
the IRAC-selected, multi-wavelength data fusion catalogue 
in the FLS (hereafter the FLS 'Data Fusion Catalogue') 



Table 1. Composition of FLS master radio catalogue. We in- 
dicate the total number of sources in the master catalogue, the 
number with cross-identifications in the FLS Data Fusion Cata- 
logue and rcdshifts, and the number of sources with redshifts and 
SPIRE/250 /urn detections 



Total number of radio sources 



1907 



with FLS Data Fusion XIDs and known redshifts 885 
with SPIRE/250 /im and known redshifts 436 



presented by Vaccari et al. (2011 in prep.). The Data Fu- 
sion Catalogue is a S'pifeer/IRAC-selected wide-area multi- 
wavelength catalog covering the ~ 60 deg 2 extragalactic 
fields covered by Spitzer/IRAG and MIPS 7— band imaging. 
The main selection of the catalog requires an IRAC 3.6 or 
4.5 /im detection, since the two Spitzer channels reach about 
the same depth. MIPS 24 /im detections are associated with 
IRAC sources to improve their positional accuracy, and the 
MIPS 70 and 160 detections are confirmed by a MIPS 
24 fim detection to increase their reliability. 

In this paper we use the version of the Data Fusion em- 
ployed in HerMES SDP wor k. For the FLS fi eld, we thus 
use the IRAC catalo g from lLacy et ail l|2005h . the MIPS 
24 /im catalogue from lFadda et al.l (|2006h , and MIPS 70 and 
160 /im catalogues produced by the HerMES tea m using the 
SSC provided software (e.g. iFraver et al.l [200 9). We com- 
bine the mid and far- infrared data from Spitzer with optical 
data (ugriz) from the Isaac Newton Telescope (Solares et 
al., 2011, in prep.) as well as redshift information from the 
literature. 

The redshifts come from the Sloan Digital Sky Survey 
(SDSS) spectroscopy and photometry as well as dedicated 
follow-up of many radio and mid /far-infrared selected tar- 
gets by several g roups ( e.g.. iMartmez-Sansigre et al.ll2005l: 
iPapovich et all 120061 : IWeedman et al.l 1200(1 " lLacv et al. 
120071; iMarleau et all 120071; lYan et all 120071; ISaiina et~ 
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120081 ; fPasvra et al.| [2009). As the photometric redshifts from 
the SDSS do not extend accurately above z — 1, higher 
redshift sources will be dominated by the selection crite- 
ria of these different groups. We can compare the optical 
magnitudes and mid-IR flux densities of the sources with 
and without known redshifts. We find that around 100 ra- 
dio sources with known redshifts are not detected in the 
z— band, but are detected at 24 /im at brighter flux densi- 
ties than most sources without redshift information. Hence, 
as faint z— band sources typically lie at higher redshifts, then 
this observation is consistent with the specific targetting of 
bright 24 /mi sources for spectroscopic follow-up at high red- 
shift. We discuss how we deal with this selection in §4. 

We cross-correlated the radio catalogue with the 
FLS Data Fusion Catalogue using a 2arcsec search ra- 
dius between the radio and mid-IR (3.6 ( im) positions. 
Exten ded /multi-component sources from ICondon et al.l 
(2003) were inspected by eye and five were reclassified as 
being two or more separate sources due to the presence 
of more than one optical/near-IR counterpart to individ- 
ual radio components. We therefore obtained a master cata- 
logue of 1907 radio sources of which 885 have spectroscopic 
or photometric redshifts from the Data Fusion Catalogue 
(see Table [1]). We illustrate in Fig. [1] the distribution in 
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(l + z) 

Figure 1. Redshift/radio luminosity distribution of 885/1907 ra- 
dio sources in our master catalogue with known redshifts. The 
red symbols within the dashed rectangles indicate our moderate 
and high redshift sub-samples represented by squares and circles 
respectively. Note the sub-samples are chosen in redshift ranges 
where they are likely to be most complete (see Fig- EJ . 

redshift/luminosity space of the sources from the master 
catalogue with known redshifts. Our search radius and the 
sky density of the FLS Data Fusion Catalogue imply that 
12/1571 (i.e. < lper cent) of our cross-identifications are by 
chance. 

While the redshift information for our sample is incom- 
plete, it is only important for sources that potentially satisfy 
our radio luminosity selection criteria and are hence included 
in our radio-loud sample. However, in the subsequent sec- 
tions we present the selection of our radio-loud AGN samples 
in two different redshift ranges, assess how complete these 
are by comparisons to models based on the known evolution 
of the high redshift radio-loud population, see §2.3, and how 
this selection will effect our sample, see §4. 



jYun fc Carillill2002l ) and therefore a S FR of ~ 6000 Mgyr" 1 
using the relations of lKennicutj l|l998h . Hence, this luminos- 
ity would be extreme for a starburst galaxy. The upper limit 
is imposed as radio sources with luminosities greater than 
this cut are rare in the volume probed in this study. We find 
one source with such a luminosity (L1.4 ~ 10 27 ' 5 WHz -1 at 
z ~ 2, see Fig. [TJ which is identified as a SDSS QSO. We 
consider it no further in this study, but note that this radio- 
loud QSO is not detected in our SPIRE observations. We 
also find that all our 'radio-loud' AGN would also be classi- 
fied as radio-loud by the rest-fr ame 5 GHz to B— band flux 
ratio according to the criteria o f lKellerman et al.l (l989) bar 
three sources in the high redshift bin which have ratios just 
below the cut-off value of ten. 

We then separate the luminous radio sources into mod- 
erate (0.4 < z < 0.9) and high (1.2 < z < 3) redshift 
samples with 15 and 16 sources respectively (out of a to- 
tal of 36 radio sources from the master catalogue with 
25 < log(Li. 4 /WHz _1 ) < 26.5). We chose these two red- 
shift bins since the redshift distribution of the luminous ra- 
dio sources peaks in these ranges (see Fig. [T]) and hence 
we should obtain the most complete sub-samples possible 
given the data available (see below for estimates of their 
completeness). We note that the general decrease in known 
redshifts at z ~ 1 seen in Fig. [T] is due to the ineffective- 
ness of SDSS photometric and spectroscopic redshift esti- 
mation above this redshift. Hence, all the sources in the 
moderate redshift sample have redshifts from SDSS (4/15 
are spectroscopic with ther remainder being photometric). 
Sources with higher known redshifts are generally from tar- 
getted follow-up of various classes of object as well as the 
occasional SDSS QSO. All the redshifts in the high redshift 
bin are spectroscopic and come from these various follow-up 
projects. Interestingly, these two redshift ranges also cover 
similar length cosmic epochs of about 3 Gyr each. The me- 
dian radio luminosities of both sub-samples are very similar: 
log(Li.4/WHz _1 ) = 24.9 and 25.0 for the moderate and high 
redshift samples respectively. 



2.2 Radio-loud Selection & Sub-samples 

To obtain accurate luminosities, radio spectral indices are 
required, so we cross-corr elated the master catalogue with 
the 610 MHz catalogue of iGarn et~ai1 (|2007T ) finding coun- 
terparts within 6 arcsec for 68 per cent of the master sam- 
ple. We use a 6 arcsec search radius to account for the po- 
sitional accuracy of the 610 MHz data. For radio sources 
without 610 MHz counterparts we assumed a spectral in- 
dex with a value of a = —0.75 (Si, oc u a ) consistent 
with the mean value found for faint radio sources in gen- 
eral (AGN and starbursts alike e.g. Ilbar et all 120091 ). We 
note that the sample here has a slightly steeper mean ra- 
dio spectral index (a = —0.82), but the relative limits 
of the 1.4 GHz and 610 MHz survey result in bias against 
sources with a flat spectrum at low flux densities. We select 
our radio-loud AGN sample with luminosity density cuts of 
25 < log(Li.4/WHz _1 ) < 26.5. The lower limit is chosen to 
ensure our sources are genuinely radio-loud and to minimise 
the number of extreme star forming galaxies (SFG) selected. 
Indeed, this lower radio luminosity is equivalent to a total IR 
(8 — 1000 fim) luminosity of ~ 3 x 10 13 Lq from the correla- 
tion of far-IR and radio luminosities for star forming galaxies 



2.3 Completeness 

In Fig.f2]we show the observed distribution of radio luminosi- 
ties in each redshift sample, and compare them to the mod- 
eled luminosity distributions over the same volume derived 
from the SKA Simulated Skies (S-cubed, IWilman et al.l 
l200ct ) at the radio flux density limit of the FLS (0.115 mjy). 
As well as the total number of sources predicted in these 
luminosity redshift bins, we also indicate the number of ex- 
treme SFGs (SFR> 6000M Q yr _1 ) predicted. The class of 
AGN from S-cubed which dominate this distribution ar e 
the low-luminosity radio-loud AGN (|Wilman et al.1 12008T ). 
The evolution of th is population is taken from 'model C of 
IWillott et all [|200lD and is reasonably well constrained up to 
z — 2. We then apply a high redshift decline in space density 
re presented by ( 1 + g) above z = 2.5 as recommended 
in lWilman et alj (|2008l ). There is also a small, ~ 6 per cent, 
contribution to the number of sources predicted by S-cubed 
of 'radio-quiet' AGN whose evolution is less well constrained 
by observation. We have included a 30 per cent uncertainty 
in the predicted number of radio sources from S-cubed to 
represent the uncertainty in the evolution of the luminosity 
function for the low luminosity radio-loud AGN population, 
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26.5 



log(L 14 /WHz ') 



Figure 2. Observed number distribution versus radio luminos- 
ity density of sources in our moderate and high redshift samples 
(black solid histograms). The dashed line represents the distribu- 
tion of the total number of radio sources expected in this volume 
from the SKA Simulated Skies (S-cubed, Wilman et al. 2008), 
where the number expected to be SFGs is indicated by the dot- 
dashed line (none are predicted in the moderate redshift sam- 
ple). The shaded region represents a 30 per cent uncertainty in 
S-cubed. In comparison to S-cubed, our moderate redshift sam- 
ple is 100 per cent complete and our high redshift sample is 14 per 
cent complete. 



in particular the high redshift cutoff and the less well con- 
strained 'radio-quiet' population, as well as sample variance 
for a survey field covering only a few square degrees. 

We find that our moderate redshift sample is complete 
given the uncertainties we ascribed to S-cubed. However, 
we find that the number of sources predicted by S-cubed 
exceeds the number we observe in the high redshift bin im- 
plying a 14 per cent completeness. The numbers of sources 
we find in each sub-sample compared to the number pre- 
dicted from S-cubed is given in Table [2] The number of 
sources deficient in our high redshift sample (and at other 
redshift ranges) can be accounted for by the lack of redshift 
information in the master sample (Table [TJ. We account 
for any bias in our samples, e.g. mid-IR selection of known 
high redshift radio- loud AGN, in §4 by considering the full 
limits of the completeness and we demonstrate that we can 
still constrain the range of mean SFRs for these samples by 
making two extreme assumptions about the sources missing 
from our sample. 




redshift=2.31 



redshift=0.64 




1 10 100 1000 

rest— frame wavelength (/zm) 

Figure 3. Example SED fits to the available IR photometry 
from 3.6 to 500 (im where we show rest-frame luminosity plotted 
against rest-frame wavelength. The red lines indicate the best fit 
starburst template and the range of templates within < 1 

and the black lines indicate the maximum normalisation of the 
AGN template to the lowest mid-IR photometry. The filled cir- 
cles indicate the Spitzer and Herschel photometry used in the 
fitting. Note that in most cases the uncertainties are smaller than 
the symbols. We present an object from both the low redshift 
sample (lower panel, z = 0.645, L IR = 5.77 ± 0.58 X 10 11 Lq) 
and the high redshift sample (upper panel, z = 2.31, L;n = 
1.84 ±0.17 x 1O 13 L ). 



3 ANALYSIS 

3.1 IR Luminosities 

We extracted SPIRE flux densities at the positions of all 
radio and 24 ;xm sourc es using the HerMES XID method 
l|Roseboom et alj|20ldh . This approach minimizes the effect 
of source blending, as the SPIRE flux densities are estimated 
via linear inversion methods using the positions of known 
24 /im sources, or radio position if there is no 24 /im coun- 
terpart, as a prior. In Roseboom et al. 2010 the 250 /im flux 
density uncertainty is estimated to be 7.45 mjy from injec- 
tion and recovery of mock sources into the observed maps. 
Flux density uncertainties are obtained from the RMS of 
input-output flux densities and consequently include con- 
tributions from both instrumental and confusion noise. We 
find 436 sources having 250 /J,m counterparts with > 3 a de- 
tections and known redshifts. For the radio-loud sources we 
find 4/15 and 9/16 with significant 250 (im detections in the 
moderate and high redshift bins respective!^. 

The relative depths of the 24 /im and 250 /im data 
available for this field have some bearing on the how the 
250 /J,m sources are found by the XID method. The 24 /im 
imaging data of this field is relatively shallow with respect 
to the SPIRE data (compared to other fields to which 



2 For reference, 24 per cent of the radio sources with unknown 
redshifts have significant detections in the SPIRE wavebands 
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this method has been applied). Hence, there may be non- 
negligible 250 fim flux remaining in the field which has not 
been extracted due to the lack of a 24 /im (or radio) coun- 
terpart. We visually inspected all 250 /im detections of the 
radio-loud sources in the SPIRE image and they all appear 
isolated with no sources close enough to them which could 
significantly effect the measurement of their SPIRE flux den- 
sity. 

We derive total (8 — 1000 pirn) IR luminosities by fit- 
ting all the data available for the 436 radio sources across 
the Spitzer /IRAC+MIPS and He rschel/SPIRE bands fol- 
lowing the method outlined in ISvmeonidis et ail (2009, 
see Fig. [3]). In all cases we use the Spitzer / '24 (j,m and 
i/ersc/ieZ/250/350/500 /im photometry although in some 
cases the 350 and 500 fjxa photometry have extremely large 
uncertainties due to their low SNR, < 3 and do not signifi- 
cantly affect the values of y 2 dervived. This fitting metho d 
uses all the models from ISiebenmorgen fc Kriigell (|2007t ) . 
which cover a wide range of SED types, and finds the best 
fit using standard \ 2 minimisation from which a total IR 
luminosity is calculated. Uncertainties in the IR luminosity 
are derived from the range of values obtained from SED fits 
which differ from the best fit by Ay 2 = (x? — Xmin) < 1- 

3.2 AGN contribution to the far-IR luminosity 

A further issue to consider, if we are to use the total IR 
luminosities as indicators of SFR, is the AGN contribution 
to this luminosity which could lead to an over-prediction of 
the star formation rates. This issue is especially important 
because o ur sources are sel e cted to be AGN. In a similar 
fashion to lSvmeonidis et all (|20101 ) , we addres s this i ssue by 
normalising a QSO template from Elvis et al. ( 1994) to the 
data point with the lowest luminosity from our photomet- 
ric dataset of 3.6 — 24 /im, as the AGN emission must be 
constrained by our photometry. If we use ot her AGN SED 
mode ls (e.g. type 1 and type 2 AGN from iPolletta et al.l 
120071 ) we find that our estimates of the upper limits to the 
AGN luminosities and ratios of AGN to total IR luminosi- 
ties change little, < 10 per cent (and therefore even less for 
the final SFR). Su ch model SEDs are broadly similar to the 
lElvis et al.l l| 19941 ) templates which in the IR are generally 
flat (in vL v ) out to the far-IR where they then drop sharply. 
We note that only fou r of our high reds hift sources have mid- 
IR spectrsocopy from lYan et alj ( |2007l ). hence we do not use 
these data to constrain SED fits. 

We then estimate the AGN contribution to the total 
infrared luminosity by integrating the QSO template in the 
8 — 1000 /im region and subtract this from our total infrared 
luminosity to obtain a star foming IR luminosity for each 
object. We can the n convert th i s star forming IR luminosity 
to a SFR using the iKennicuttl (|l998T ) relation. In Fig. 0]we 
show the AGN IR luminosity and ratio of AGN to total 
IR luminosity as a function of radio luminosity density for 
the radio sources detected by SPIRE in our two redshift 
samples. The IR AGN luminosity has a large scatter which is 
largely due to the moderate redshift sub-sample having lower 
IR AGN luminosities (~ 10 n Le) than the high redshift 
sub-sample (~ 10 12 £©), although we observe no trend with 
radio luminosity within a sub-sample. 

We suspect that the greater AGN IR luminosities of the 
sources in the high redshift sample is most likely due to a 
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log(L 14 /WHz l ) 

Figure 4. Upper limits to the AGN IR luminosity estimated from 
the normalisation of an AGN SED to the lowest mid-IR luminos- 
tiy, and the resulting upper limit to the ratio of AGN to total IR 
luminosity, both plotted as a function of radio luminosity den- 
sity. These sources are from both our moderate and high redshift 
samples (open squares and circles respectively). 

bias in the redshift identification toward sources with bright 
24 /im flux densities (^ 1 mjy) as discussed in §2.1. Addi- 
tionally, the flux limited nature of the Spitzer and Herschel 
data mean that the SPIRE observations are more sensitive 
to lower IR luminosities at lower redshifts. As these inden- 
tified sources comprise just 14 per cent of the high redshift 
sample they are not likely to be representative in terms of 
their AGN fraction. 

The ratio of AGN to total IR luminosity tends to be low, 
under 0.3 bar one source, co nsistent with the results seen in 
lHatziminaoglou et alj (|201Ci) . and averages around 0.15. As 
a check we apply the si multaneous AGN/starburst te mplate 
fitting routine used by lHatziminaoglou et all (|2010h to the 
radio-loud AGN studied here and we find similar total IR 
luminosities and AGN fractions. Therefore, the final SFRs 
we derive are not very sensitive to our choice of model star- 
burst and AGN SEDs. Our assumption that the mid-IR is 
completely dominated by the AGN, while conservative, also 
does not have a strong effect on the final SFR due to the 
low AGN fraction. 

3.3 Comparison between radio and IR 
luminosities 

We calculate the total IR luminosities of all SPIRE detected 
sources in order to confirm our method of measuring these 
luminosities by comparison with the radio/far-IR correla- 
tion seen in local s tar forming galaxies and now confirme d 
at higher redshifts (|Sevmour et al.ll2009l : [ivison et alJl201Ch . 
Additionally, by extrapolating this empirical correlation to 
higher luminosities we can assess the contribution of star 
formation to brighter radio sources. Our radio luminosity 
selection would be equivalent to a SFR of ~ 6000 MQyr -1 
for a pure SFG, but potentially there could be a few sources 
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Table 2. Composition of radio-loud AGN sub-samples. For both 
redshift sub-samples we present the number of sources pre- 
dicted from S-cubed, the total number found, the number with 
SPIRE/250 fim detections, the mean 250 fim flux densities of the 
detected and undetected (via stacking techniques) sources, the 
mean SFR of the detected sources, the inferred mean SFR of the 
undetected sources (assuming it scales directly with the mean 
250 fim flux density), the total mean SFR of all observed sources 
and the range of mean SFRs given the number of sources pre- 
dicted by S-cubed. 



description 




sub-sample 
moderate high 


S-cubed number predicted 

total number found 

with SPIRE/250 fim a > 3 


16 
15 

4 


116 
16 

9 


/ odctcctcd\ 
\°250 / 
/ cundctcctcd \ 
\°250 / 


(mjy) 
(mjy) 


27.2 ±2.5 
2.0 ±0.8 


39.2 ± 2.5 
6.5 ± 1.0 


( SFR 250° CtOd ) 

<SFR§o al ) 


(Meyr- 1 ) 
(Meyr- 1 ) 
(M yr-i) 


92 ±28 
6.7±2.8 
29.5 ± 11.6 


914 ± 274 

153±23 
581 ± 143 


range of (SFR) 


(Meyr" 1 ) 


18 - 41 


80 - 581 



log(L TR /Lj 



Figure 5. Total infrared luminosity plotted against radio lumi- 
nosity for all the FLS radio sources with redshifts and a > 3 a 
detection at 250 fim (including confusion noise). The sample of lu- 
minous radio sources, L1.4 > 10 25 WHz -1 , used in this work are 
indicated by open symbols and the less luminous radio sources by 
asterisks. Note the IR luminosities of the luminous radio sources 
do not include the AGN contribution (see §3.3). We fit the ob- 
served correlation in luminosities for £1.4(3 .f/z < 10 24 WHz -1 and 
Lir < 10 12 5 Lq and derive q IR = 2.40 ±0.19. Two radio-loud 
AGN lie above, but within 2 a of this correlation, however they 
have radio spectral indices inconsitent with star formation. Hence 
we conclude that their radio emission is dominated by AGN pro- 
cesses. 



with higher SFRs within the volume probed here (see S- 
cubed predictions in Fig. [2J. 

We find a strong correlation between the radio and 
IR luminosities, particularly below Li. 4 = 10 24 WHz -1 and 
Lyk — 10 12 ' 5 Lq (see Fig. [5]) which we use to verify our IR 
luminosities. Note the IR luminosities of the radio sources 
in the moderate and high redshift samples have had the 
AGN contribution removed. We define the ratio of radio to 
IR luminosity as giR = log(LiR,/ X1.4) + 14.03 (as used in 
ISaiina et alfeoOSl . this definition is an equivalent, but more 
convenient form than the classical one of Helou et al., 1985). 
By fitting the correlation over these luminosity ranges we 
get a value of qm = 2.40 ± 0.19 using a biweight estima- 
tor (|Beers et al] 1990 | ). in pood ag reement with the value 
found locally (lYun fc Carillil [20021 ) and at higher redshift 
(|lvison et alJl201Ch . 

We then assume that this relation holds to higher lumi- 
nosities (i.e. to SFRs > 6000 M^yr -1 ) and note that 2 radio 
sources with log(Li.4/WHz _1 ) > 25 lie just within 2 a of qm 
The proximity of these 2 sources to the correlation may mean 
that these sources have a non- negligible contribution of star 



formation to their radio luminosity. The AGN fraction of the 
total IR luminosities for these sources is low, ^ 10 per cent. 
However, upon closer inspection these two sources have ra- 



dio spectra which are either too steep, a; i4 



-1.78, or too 



flat, qj.4 = —0.22 , compared to the canonical value for star 
forming galaxies l|Condonlll992l ). The fraction of the radio 
luminosity due to star formation, assuming the star forming 
component lies precisely on the correlation, is 25 — 30 per 
cent. Hence, we conclude that their radio emission is dom- 
inated by AGN processes and retain them within our high 
redshift sample. 



3.4 Stacking the non-detections at 250 /im 

We can obtain an approximate constraint on the far-IR lumi- 
nosity of the radio-loud AGN not detected at 250 /im in each 
sample by employing stacking techniques to obtain mean 
250 fim flux densities for those sources. By assuming the 
same distribution of redshifts, IR SED types and ratios of 
AGN to total IR luminosity, we can argue that the mean 
SFRs of the undetected and detected samples scale directly 
with the 250 fim flux densities within both redshift ranges. 
Therefore we stacked the 11 and 7 sources not detected at 
250 fim in each sub-sample and find the mean flux densities 
reported in Table [2] The uncertainties in flux densities of 
the stacked sources are simply those of the mean. 



3.5 SFRs in local (z < 0.1) radio-loud AGN 

In order to examine any evolution of the mean SFR of radio- 
loud AGN over cosmic time we need a local baseline to com- 
pare with. Recently published Spit zer/MIPS observati ons of 
the local (z < 0.1) 3CRR sample l|Dicken et al.ll2010h pro- 
vide an excellent opportunity to assess star formation in 
the nearby radio-loud population. The 1.4 GHz luminosity 
densities of this sample, derived from the 5 GHz values in 
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log(l+z) 

Figure 6. Range of mean SFRs plotted as a function of redshift for radio-loud AGN with 25 sg log(Z,i. 4 /WHz _1 ) ^ 26.5 (shaded 
regions). At 0<z<0.1 the data is from our 3CRR local reference sample and at 0.4<z<0.9, and 1.2<z<3.0 from our moderate and high 
redshift sub-samples respectively. The open rectangles indicate the results from Hardcastle ct al., (2010) using Herschel observations of 
sources with a similar range of radio luminosities. The points with error-bars present the approximate mean SFRs of X-ray selected AGN 
over the range of rcdshifts indicated from Lutz et al. (2010, open circle) and Shao et al. (2010, asterisk). 



iDicken et all l|2010h assuming a = —0.75 and oc v a , fall 
within the 25 < log(Za.4/WHz ) < 26.5 range of our sub- 
sample selection. The 3CRR sources were selected to only 
include sources with Fanaroff-Riley class II morpholgies (i.e. 
those with radio lobes which are brightest at their edges, 
iFanaroff fc Rilevll 19741 ). However the lower radio luminosity 
density limit used in our work very closely corresponds to 
the luminosity density, log(Li.4/WHz _1 ) = 25.1, at which 
the radio-loud population switches from mostly containing 
class I sources to mostly containing class II sources. Further- 
more, this local sample is not sensitive to the low end of our 
radio luminosity density range at z — 1, and therefore may 
not be 100 per cent complete. Dicken et al. (2010) derive 
rest-frame 70 fim luminosities from their Spiteer/MIPS ob- 
servations which they compare with the [OIII] emission line 
luminosities of the local 3CRR sample. They find a broad 
correlation implying that generally the 70 /im luminosity is 
due to the AGN. However some 3CRR sources, which show 
evidence of star formation from their optical spectra, gen- 
erally lie above this correlation, i.e. they have an excess of 
70 /im luminosity compared the [OIII] emission. These au- 
thors postulate that this 70 fim excess could be due to star 
formation. 

Here, we estimate the range of mean SFR in this sample 
using two assumptions. To obtain an upper limit we assume 
that all of the 70 /jm luminosity is due to star formation. 
To obtain a l ower limit we use the linear regression fit by 
IDicken et all (|2010T) to the correlation of the OIII and 70 /im 
luminosities to estimate the AGN only 70 /im luminosity. We 
then subtract the AGN luminosity from the total 70 /im lu- 
minosity for all sources lying more than 0.3 dex above the 
correlation in order to obtain a starburst only 70 /im lumi- 
nosity. In both cases we convert the 70 fj,m luminosities to 



total I R luminosities using the relation of ISvmeonidis et al.l 
( 2008) and then to SFRs using the Kennicutt (1998) relation 
as before. Due to the size of the sample and the influence of 
one very luminous source we use the median inferred SFR 
and find that the range of typical SFRs for the local 3CRR 
sample is 3.4 — 4.2M0yr _1 from these two assumptions. 



4 RESULTS 

In Table [5] we report the mean SFR, (SFR) , of the radio- 
loud AGN detected at 250 /im in each sub-sample. The SFRs 
of individual sources are derived from the total IR lumi- 
nosities, minus the A GN contrib u tion ( see § I3.2jl , using the 
conversion factors of iKennicutt] |l998) . We find values of 
92 ± 28 M yr _1 and 914 ± 274 M©yr _1 in the moderate and 
high redshift bins respectively. For the sources undetected 
at 250 /im we find stacked 250 /im flux densities which are a 
factor eleven and seven lower than the mean flux densities of 
the detected sources (see Table [5| for the moderate and high 
redshift sub-samples respectively. It is unsurprising that un- 
detected sources have a mean flux density lower than those 
detected, but the fact they are considerably lower (i.e. not 
just below our 3 a cut) suggests that these radio-loud AGN 
have a wide range of intrinsic SFRs. We report the SFRs 
of the undetected sources in Table [2] obtained from the ra- 
tio of the mean 250 /im flux densities of the detected and 
undetected sources and the measured SFR of the detected 
sources. Then we estimate the total mean SFR in each subset 
by combining the mean SFR of the detected and undetected 
sources weighted by the number in each group. The esti- 
mated total mean SFRs for the total sample are therefore 
29.5 ± 11.6 M Q yr _1 and 581 ± 143 Mgyr" 1 for the moderate 
and high redshift bins respectively. 
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The moderate redshift sample is complete within the 
uncertainties of the S-cubed simulation (we find 15/16 pre- 
dicted sources in this redshift/luminosity density paramter 
space). Hence, we can directly calculate the mean SFR of 
the low redshift sample by summing the observed SFRs and 
dividing by the number sources. The uncertainties are sim- 
ply those of the measured SFRs, which directly come from 
the uncertainties in the IR luminosities, combined with the 
30 per cent uncertainty in the S-cubed model. As the lat- 
ter are so much greater than the former our uncertainties 
are dominated by the conservative uncertainties we used 
in S-cubed. We find a mean SFR for this sub-sample of 
29. 5± 11.6 Moyr -1 which is equivalent to the range of values 
of 18 -41 M yr _1 . 

As we saw from Fig. [2] the high redshift sub-sample 
is incomplete, although we can quantify the incompleteness 
from comparisons to the S-cubed simulation. The number of 
radio- loud AGN expected from S-cubed is given in Table [2] 
We cannot estimate the properties of sources not included in 
our high redsift sample due to lack of redshift information. 
However, we can estimate likely lower and upper limits on 
the mean SFR from two simple assumptions. Firstly, to es- 
timate the lower limit we assume that all the sources missed 
have SFRs of zero and then scale the mean SFR by the in- 
completeness (i.e. the lower limit is x the mean SFR for 
the observed fraction). While 24 per cent of the sources with 
unknown redshifts have 250 /jm detections we have no way 
of knowing how many of these fall into our high redshift sub- 
sample, hence this method of determining our lower limit is 
the most robust approach. Secondly, for the upper limit we 
assume that all sources not included have mean SFRs iden- 
tical to the detected fraction, i.e. the upper limit is simply 
the measured mean SFR for the detected fraction. Hence, 
we calculate the range of mean SFRs for the high redshift 
sub-sample to be 80 — 581 Moyr - . 

We compare these constraints with those found for the 
local 3CRR sample and the recent results of Hardcastle et al. 
(2010) in Fig.|6]who measure IR luminosities from Herschel- 
ATLAS observations of radio sources occupying a similar 
region of redshift/luminosity parameter space. We see an 
increase in the mean SFR of radio-loud AGN with cosmic 
look back time. In the local Universe we found the mean SFR 
of z < 0.1 radio-loud AGN to be 3.4 - 4.2M yr -1 , whereas 
at moderate redshifts, 0.4 < z < 0.9, we constrain it to be 
~ 5 — 10 times greater and in our high redshift sample we 
find it to be ~ 20— 150 times greater. While these ranges of 
mean SFRs are wide we observe a clear trend of increasing 
mean star formation rate with redshift in radio-loud AGN in 
the luminosity density range 25 < log(Z/i.4/WHz -1 ) < 26.5, 
a trend t hat is also seen over a sm aller redshift range in the 
results of lHardcastle et alj l|2010l ). 

We can quantify this rate of increase by fitting a 
straight line through the shaded regions of Fig. [6] via lin- 
ear regression. We then find that the mean SFR of radio- 
loud AGN in this luminosity range evolves as (1 + z)®, 
where we measure the value of Q = 4.2 ± 0.8. This 
value for the evolution is strong and greater than that 
measured for the evolution of the star forming luminosity 
function (which typically has values o f Q ~ 3 as traced 
by IR surveys iLe Floc'h et al l 120051; iHuvnh et ail 120071 ; 
iMagnelli etafl 120091 ; iRodighiero et all l201Ch . We can also 
compare our results with the mean SFRs of high redshift 
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Figure 7. Mean star formation rate, (SFR), of the host galaxies 
of radio-loud AGN plotted as a function of radio luminosity den- 
sity for each of the redshift sub-samples. The shaded regions for 
each bin then represent the range of range of (SFR) assuming ei- 
ther (i) all the sources missed have (SFR) = (the lower limit) or 
(ii) all the sources missed have (SFR) equal to the sources found. 
The dashed lines represent the range of (SFR) for the whole of 
each redshift sample as given in table [2] 



AGN selected at other wavelengths. The mean SFRs of X- 
ray selected A GN, L2-iokeV > 1 43 e rgs -1 , have been stud- 
ied recently bv lShao et al.l 1 2010T ) and lLutz et al.l (2010) who 
find that such sources have mean SFRs within, but at the 
low-end of, the range of values found our high redshift bin. 
We illustrate those results in Fig. [6] using the same Kenni- 
cutt total IR luminosity to SFR conversion as before and 
converting the Shao et al. 60 /im monochromatic luminos i- 
ties u sing the formula presen t ed inlSvmeonidis et al.l (2008). 
Also, lHatziminaoglou et al] (|2010f ) find a similar range of 
SFRs for a heterogenous sample of AGN above 2 = 1, sug- 
gesting that this increase is common to different types of 
AGN activity. 

If we sum the observed star formation in each red- 
shift bin we can calculate the comoving star formation 
rate density due to the host galaxies of the radio-loud 
AGN in each redshift sub-sample. We find values of ~ 
2.5 x 10 -5 Moyr -1 Mpc -3 for the moderate redshift bin and 
1 - 5 x 10 -4 M yr -1 Mpc -3 for the high redshift bin. For the 
local redshift bin, the star formation density due to the host 
galaxies of the radio-loud AGN is ~ 4x 10 -8 M Q yr -1 Mpc -3 . 
We can compare these SFR densities with the globally mea- 
sured SFR history from a variety of different methods (e.g. 
Hopk ms fc Beacomll2006h . We observe that the relative con- 
tribution of the host galaxies of radio-loud AGN to the 
total comoving SFR density increases with redshift from 
~ 0.0004 per cent in the local sample to ~ 0.03 per cent 
and ~ 0.1 — 0.5 per cent for the moderate and high redshift 
samples respectively. 

In Fig. [7]we show the (SFR) as a function of radio lumi- 
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nosity density for each of our two redshift sub-samples. We 
calculate upper and lower limits for each luminosity density 
bin as we did for the whole sample. The upper and lower 
limits are indicated by the grey shaded regions. Note due 
to the fact we only detect 4/15 sources in the moderate 
redshift sample we have to increase the bin size by a fac- 
tor of three compared to the high redshift sample. We also 
overlay the upper and lower limits for the whole of each sub- 
sample as indicated by the dashed lines. We see no evidence 
for any trend of mean SFR with radio luminosity for either 
sub-sample, although the constraints for the highest radio 
luminosity density bin of the high redshift sample are not 
so strong. 



5 DISCUSSION 

We observe that radio-loud AGN in the distant Universe 
have an increasing mean SFR with cosmological look back 
time in the 25 < log(Li.4/WHz _1 ) < 26.5 radio luminos- 
ity density range. In the local Universe, z < 0.1, the mean 
SFR of the 3CRR sample is 5 — 10 times less than that in a 
moderate redshift sample, 0.4 < z < 0.9. We note that the 
3CRR sample was also selected on FRII radio morphology 
which suggests we may not be comparing identical popula- 
tions, and it may not be 100 per cent complete. Another re- 
cent study has examined the IR luminosities of bright radio 
sources with Herschel- ATLAS observations of the GAMA-9h 
field. With a similar radio luminosity cut as our moderate 
redshift sub-sample, Hardcastle et al. (2010) find a mean 
SFR of between 20 and 5OM0yr _1 , increasing across our 
moderate redshift bin (see Fig.[S|. This range of mean SFRs 
is consistent with that found here, 18— 41 Moyr -1 , allow- 
ing for the slightly different source selection, the different 
method of estimating IR luminosities and the fact these au- 
thors do not subtract any AGN contribution to the total IR 
luminosity. 

We find the increase in mean SFR of radio-loud AGN 
hosts (parameterised as ~ (1 + z)®, Q = 4.2 ± 0.8) to 
be greater than that of the IR luminosity function which 
traces the evolution of the general star forming popula- 
tion. This greater rate of increase with redshift, compared 
to the regular star forming population, suggests that some 
of the star formation may be directly associated with the 
radio-loud AGN activity. The increase of mean SFR with 
redshift of AGN is also seen in X -ray selected AGN (e.g. 



iLutz et "all l20ld; | Shao et all 20101) and in a heterogenous 
sample of AGN l|Hatziminaoglou et alj |2010). Alternatively, 
our results could reflect an increase in the stellar mass of 
the host galaxy, since high stellar mass galaxies have SFRs 
which increase with redshift (e.g. lJuneau et al.ll2005h. This 
interp retation would fit in with the recent Tadhunter et al.l 
(j2O10T ) result who find that at low redshifts, z < 0.7, not 
all ULIRGs are massive enough to host radio-loud AGN. 
If the stellar masses of ULIRGs increase with redshift then 
ULIRGs would be more likely to host radio-loud AGN at 
higher redshifts. 

While it is likely that the redshift information for the 
high redshift sample is biased toward sources that have 
bright 24 /im flux densities (see § 12. ip . our approach of deter- 
mining a range of mean SFRs given two extreme assump- 
tions alleviates much of the concern about selection bias. 



The remaining principle source of uncertainty is the S-cubed 
model, used to quantify how complete our sub-samples were. 
As discussed earlier, our uncertainties in S-cubed are very 
conservative. S-cubed treats the AGN and SFGs as separate 
populations, i.e. it does not include hybrid radio sources 
exhibiting both processes simultaneously. We can thereby 
compare the expected number of radio-loud AGN regard- 
less of whether there is ongoing star formation in their hosts 
or not. The contribution to the SFR density of the host 
galaxies of radio-loud AGN in the high redshift bin is in- 
teresting as the SFR density at this epoch i s dominated by 
LIRG s and ULIRGs |Le Floc'h et al.ll2005l ; ISevmour et"atl 
2010). As 0.1 — 0.5 per cent of the SFR density consists of 
LIRGs and ULIRGs which host the radio-loud AGN, we can 
infer a duty cycle of 0.001 — 0.005 for radio-loud AGN activ- 
ity in such sources, assuming that each LIRG and ULIRG 
goes through at least one radio-loud phase. The typical time- 
scale of a radio- loud phase of an AGN is around ~ 10 Myr 
for extended radio sources <|Milevlll980h and likely shorter 
for the less luminous sources with smaller radio lobes consid- 
ered here. Given this lifetime and the estimated duty cycle 
of 0.001 - 0.005 we can estimate that LIRGs and ULIRGs 
undergo a radio-loud AGN phase every 2 — 10 Gyr. Hence, 
during the 3 Gyr time span covered by the high redshift sub- 
sample we could expect perhaps one major phase of radio- 
loud AGN acti vity at a rate similar to that expected from 
major mergers (|Hopkins et al.ll2010T ). 

The feedback models which que nch star formation 
by evoking a ra dio- loud phase (e.g. ICroton et alj 120061 ; 
iBower et al.ll2006h are most important at late times, i.e. be- 
low z < 1, but they must occur at higher redshifts in order 
to prevent the most massive galaxies, formed at early times, 
from growing significantly more. However, in this work we 
observe many AGN in our high redshift sub-sample which 
are in a state equ i valent to the 'radio -mode' feedback of 
ICroton et~ al. (2006); Bo wer et al] |2006h and simultaneously 
have very high SFRs while feedback processes are predicted 
to be occuring. 



6 CONCLUSIONS 

We have examined the incidence of far-IR emission and 
inferred SFR of luminous radio-loud AGN in a moderate 
redshift, 0.4 < z < 0.9, and a high redshift sub-sample, 
1.2 < z < 3, as well as a local, z < 0.1, comparison sample. 
We have: 

• constrained the mean SFR of radio-loud AGN to be 
3.4-4.2, 18-41 and 80-581 M yr _1 for the local, moderate 
and high redshift samples respectively, hence, we measure 
the evolution of the mean SFR to be ~ (1 + z) 4 ' 2±0 ' 8 , 

• observed no strong trends of SFR with radio luminosity 
in any redshift bin, 

• estimated that the host galaxies of radio-loud AGN in 
the high redshift sub-sample contribute 0.1 — 0.5 per cent to 
the total SFR density at that epoch and if all LIRGs and 
ULIRGs have a radio-loud phase we infer a duty cycle of 
0.001 — 0.005 in such sources. 

These results demonstrate that in the distant Universe 
a considerable amount of star formation is occuring in galax- 
ies hosting a radio-loud AGN, consistent with the frequent 
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evidence for high SFRs in classic high redshift radio galaxies. 
The mean SFR evolves more quickly than the IR luminosity 
function implying that some of the star formation is directly 
related to the radio-loud AGN activity. Both starburst and 
active nuclear processes have relatively short time-scales so 
their co-existence in many objects suggests that bursts of 
star formation and jet activity either are quite common or 
connected via 'feedback'. But is the jet initiating or quench- 
ing star formation, or are the processes independent? We 
cannot answer such questions here, but we shall be able to 
do so with follow-up of individual sources (to search for out- 
flows of jet-triggered star formation or for mergers triggering 
both) and with the huge sample that will be provided by the 
full HerMES data set combined with improved redshift in- 
formation. 
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